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ON SOME INTEGRALS INVOLVING THE HURWITZ ZETA 
FUNCTION: PART 2 



OLIVIER ESPINOSA AND VICTOR H. MOLL 



Abstract. We establish a scries of indefinite integral formulae involving the 
Hurwitz zeta function and other elementary and special functions related to 
it, such as the Bernoulli polynomials, lnsin(7rg), Inr(g) and the polygamma 
functions. Many of the results are most conveniently formulated in terms of a 
family of functions Afc(g) := k(^'{l — k, q), A; 6 N, and a family of polygamma 
functions of negative order, whose properties we study in some detail. 



u 

"ti I The Hurwitz zeta function, defined by 



1. Introduction 



oo _. 

(1-1) C(2,9) = Vt r 



n=0 



^ ■ for z G C, Rez > 1 and q ^ 0, —1, —2, • • • , admits the integral representation 

OO- (1-2) az,q) = — -/ - — -e-'dt, 

o 



where T{z) is Euler's gamma function, which is valid for Re z > 1 and Re g > 0, and 
can be used to prove that C(^)'?) has an analytic extension to the whole complex 
(^ • plane except for a simple pole at z — 1. 

For Ke z < 0, ({z^q) admits the following Fourier representation, originally de- 
rived by Hurwitz, in the range < q < 1: 



, X , ^ 2r(i-z) 

(1.3) C(^,g)- ^ ' 



cos{27Tqn) /""-^A V^ sin(27rgn) 



/7rz\ v-^ costZTrqnj (t^z\ v-^ 



n=l 



•i-H . 

^^ ' This representation was used in ||7| to obtain several definite integral formulae 

H , involving C(z, q). A derivation of (|l^ can be found in |1J], page 268. An alternative 

proof, based upon the representation 

n^-^ rr^ 1"°° H\ — i 

(1-4) C(^,9) = ^ + V-W T^^^dr^^. 



z-l 2 7o {t + qY+^ 

where {t\ is the fractional part of t, has been given by Berndt B. The expression 



(1.4) is employed in M to give short proofs of several classical formulae, including 



Lerch's beautiful expression 

(1.5) inr(g) = c'(o,<z)-C'(o). 
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2 OLIVIER ESPINOSA AND VICTOR H. MOLL 

In this paper we continue the work, initiated in M, on the exphcit evaluation of 
integrals involving C(z, q). Special cases of C(z, q) include the Bernoulli polynomials, 

(1.6) B^{q) = -mC(l-m,(7), m G N, 

defined by their generating function 

(1-7) ^ = T.B'rM)- 

■m— 

and given explicitly in terms of the Bernoulli numbers B\^ by 



(1.8) B„,{q) = EiJ^fe?" 

fc=0 

the digamma function. 



• -k. 



^ -C(^,g) 



(1.9) V(<?):=:^lnr(g)=lim 

dq z^i \_z — 1 

and the polygamma functions, 

(1.10) V^"^(9) = (-l)"^^"i!C(m + l,g), m e N, 
defined by 

(1.11) V'^^n?) := rf^^('^)' "^^N- 

An important property of the Hurwitz zeta function, which will be essential for 
the indefinite integral evaluations presented in Section g, is the following: 

(1.12) —az,q)=-zaz + l,q). 

The rest of this paper is organized as follows. In Section || we consider the 
evaluation of indefinite integrals of functions of the form f{q)C.iz, a + bq), using a 
simple integration by parts approach. In Section ^ we introduce and study some 
of the properties of two families of functions related to the first derivative with 
respect to the argument z of the Hurwitz zeta function (^{z, q), evaluated at z equal 
to nonpositive integers. These functions appear in connection to the indefinite 
integrals involving polygamma and negapolygamma functions, as well as In r(g) and 
lnsin(7r(7), considered in Section m These, in turn, are derived from the formulae 
obtained in Section either by direct differentiation or by taking the appropriate 
limits. Finally, in Section H we use some of the indefinite integral formulae to 
rederive some of the definite integral evaluations obtained in reference [Q and to 
present some new analogous formulae. 

2. The evaluation of indefinite integrals 

In this section we discuss a method to evaluate primitives of functions of the form 
f{q)(^{z,a + bq). This is illustrated in the cases where / is a polynomial and an 
exponential function. The resulting evaluations can be taken as a starting point to 
derive similar formulae involving other special functions in place of C(z, a + bq) . For 
instance, differentiation with respect to the parameter z leads, in view of Lerch's 
result (|l.5|), to the evaluation of primitives involving the weight lnr(a + bq), and 
thus also In sin irq, by virtue of the reflection formula for the gamma function. Also, 
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the limit z ^ m (z N leads to the evaluation of primitives involving the polyganima 
function ■0'^™~^)(g). We shall present these results, from a slightly more general 
point of view, in Section H. 

Theorem 2.1. Let ?■ e N, / be r-times differentiable and a. 6 G M. Then 

(2-1) / f{qK{z,a + bq)dq = }_^{-l) ^ ^^ ,^ _ . 



fc=i 



f^'>{q)az~r,a + bq)dq. 



b'-il-zh 
Proof. Observe that 

(2.2) ^(^^z-l,a + bq) = bil-z)C{z,a + bq), 

so that integration by parts yields 

fiq)az-l,a + bq) 



f{q)C{z,a + bq)dq = 



6(1 -z) 
1 



b{l - z) 



f'{q)C,{z-l.a + bq)dq. 



The expression (2.1) follows by repeating this procedure. D 



We now produce the evaluation of certain indefinite integrals by choosing appro- 



priate functions / in Theorem 2.1 



Example 2.2. Let n G Nq and a, 6 G R. Then the moments of C(z, q) are given by 



(2.3) / g"C(^, a + bq) dq = nlf^ ,,+i./ ^^'S ' TyC(^ - J - 1, a + bq). 

Proof. The case n = is simply the known result 

C(z-l,a + &(z) 



(2.4) j C{z,a + bq)dq 5(1 - z) " 

For n > 1, the function J{q) — q" satisfies f'^^^^\q) = n\q"^''+^/{n — k + 1)!, for 
k < n. Then ( p.l| ), with r ~ n, yields 

A (-l)fc+in!g"~fe+i 
^^ (n-fc+ 1)16^(1- z)fc ^ 

(-l)"n! 



q^Xiz,a + bq)dq ^ g (^IfeV 1),,^ - z). ^^^ ^' " + ^'^^ 

C(z — n,a + bq) dq, 



b"{l- z)r 



so that (KJ) follows from (2.4). D 



In a similar fashion we obtain: 

Example 2.3. Let m G No and a, b, c, d e M. Then 

(2.5) 



/ B„,{c + dg)C(z, a + bq) dq ^ ml J^ TtSt^t'^T^^^^'' ~J~l,a + bq). 
J 'H o^+(l- 2^)7+1 (™-j)' 
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Proof. Same as the proof for Example 3.2 with 



dq — B^i^ + dl) - ^^:^^^y B^-,+l{c + dq). 



D 



Definition. The family of fmictions ^ := {fj{q) : j G N} is said to be closed 
under primitives if for each j the primitive of fj (q) can be written as a finite linear 
combination of the elements of ^. Naturally, the family ^ is allowed to depend on 
a finite number of parameters, as in the next example. 



Example 2.4. Example 2.2 shows that 

■Sa.b '■= {Pj {Q)Ciz — m,a + bq) : j/, to G N and Pj is a polynomial in q 

of degree j } 
is closed under primitives. This follows from 



/n 
g"C(z - TO,a + bq)dq = ?i! ^ 

which is a variation of (^.3|). 



{-lyq^'^Ciz -m-l-j,a + bq) 



Example 2.5. The moments of the Bernoulli polynomials are given by 

(2.7) 

/ q-B,Aa+bq) dq ^ -^^ t ^^ (^" ^ " ^ ') B^^.^.^a^bq). 



i=o 



Proof. Use the identity (1.6) in (2.3). D 



Example 2.6. Let n G N be odd. Then 

(2.8) 

/ C(z - "■, <l)C,{z, q)dq = 2 51 /-,_ \^^ C{z ~ k, q)C{z ~n + k~l,q). 

Proof. The Hurwitz zeta function satisfies 

(2.9) g^Ci^ -n,q)^ (-1)'"'(^ - ri)fc-iC(2 - n + fc - l,q), 



SO the result follows from (2.1) with r = n, since in that case the integral on 
the right-hand side equals the one on the left-hand side, except for the prefactor 

(z-n)„/(l-z)„ = (-l)". 

D 

Note. In view of the identity 

(z-n)fc-i ^ ^^^1 {z~n)n-k 

(l-2)fc (l-z)„_fe+i' 
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it is easily seen that the terms in the sum on the right-hand side of (2.9) are equal 
in pairs, except for the central term k = r, where r G N is defined by n = 2r — 1. 
Therefore we have the alternative formula 

(2.10) J ({z - 2r + 1, g)C(z, q)dq = ^^ 2 (T- J)!'"' ^'^^ " ''' '^^ 



r — 1 

^ (i_^!:±ikzic(^ - k, q)C{z - (2r ~ fc), g). 
k=l 



(1-^)* 



Note. We have been unable to evaluate the integral in Example |2.6| for the case n 
even. Thus the question of whether the family 

(2.11) dz ■■= {(:{z-n,q)C{z~m,q): n,meN} 

is closed under primitives remains to be decided. 



Example 2.7. Let a,beM.. Then 



(2.12) e''C{z,a + bq)dq = e' 






i^iy 



3=0 



(1 - Z)j+1 



C(z- l-j,a + bq). 



Proof. Divide (2.3) by n\ and then sum over n to produce 



j eK{z.a + bq)dq = ^ ^ ^^ 



{l-z)j+i{n- i)\ 



C,{z- 1- j,a + bq). 



n=0 j=0 

The result follows by interchanging the order of summation. 



D 



Note. We have been unable to produce a finite expression for the integral in (2.12) 



Example 2.8. Let m G N. Then 



(-IP' 



(2.13) I eW^{a + bq)dq = m!e«(-l)'" ^ ^^6"--'B,(a + 6g). 

3=0 ^' 



Proof. Use the identity (1.6) and (to)j+i = (rn + j)!/(m— 1)! in (2.12) to produce 



(2.14) l'eW,-,{a + bq)dq = m!e«(-l)™+i ^ ^^6"~^Bj(a + 6g) 



j=m+ 



7 
1 ■' 



The generating function (1.7) is now employed to see that the sum from j = to 
infinity is independent of g, so it is absorbed into the implicit constant of integration. 

D 
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3. The function Ak{q) and negapolygammas 
In this section we consider the function 

(3.1) A,(g):=fc|-C(^,g) 

for fc G N. The function ^i(^) has a simple expHcit form, 

(3.2) Ai(g)=C'(0,g)=hir(g)+C'(0), 



in view of Lerch's result ( |1.5D . These functions appear in all of the formulae for 
the indefinite integrals involving the loggamma and the logsine functions studied 
in Section 0. The derivative of the Hurwitz zeta function has appeared before in 
connection to integrals of Inr(q) [g|, and in a number of related contexts, such 
as the studies of polygamma functions of negative order Q, the Barnes function 
g and the multiple gamma function iQ, and other unrelated ones such as the 
evaluation of sums of the type X]m>2 ■^™C('7i, a) Q. Here, we will rediscover, 
from a more general point of view, the intimate relationship existing between the 
functions Ak{q) and polygamma functions of negative order, first pointed out in 
reference Q. 

The nonelementary behavior of Ak (q) can be considered to be contained in the 
range < q < 1, since for q > 1 the value of Ak{q) can be obtained by repeated 
use of the following result: 

Lemma 3.1. The function Ak{q) satisfies 

(3.3) Akiq + 1) = Akiq) + kq''-^\nq. 
Proof. Differentiate both sides of the identity 

(3.4) C{^.,q)^^+az,q+l) 

q 

with respect to z at z = 1 — k. D 

As an immediate consequence of property (3.3) and definition ( ^.l| ) we have 
Lemma 3.2. For fc > 2, 

(3.5) Ak{0) = Ak{l) = kC{l-k). 

Proof. Take the limit g ^ in (|3.3[ ) and use the fact that C(z, 1) = C(-2)- D 

Lemma 3.3. For fc e N, 

(3.6) / Ak{q)dq = 0. 

JQ 

Proof This is a direct consequence of 

(3.7) / az,q)dq = 0, 

Jo 

valid for Re z < 1 . D 



HURWITZ ZETA FUNCTION 

Lemma 3.4. For fc e N, the functions Ak{q) satisfy 
(3.8) A',^,{q) = (fc + 1) 

Proof. We have 



Akiq) + ^Bkiq) 



A'k+M) - (fc + i)l:l:C(^,9) 



dq dz 



-(fc + l)|-|-C(2,9) 
-fc oz oq 



= -k 



dz 



(fc + 1) 



1 



Akiq) + -Bkiq) 



The foUowing results can sometimes be used to simphfy a formula: 
(3.9) <'(-2n) = 



-D- '^-rg^r" , ».N, 



2 (27r)2 



(3.10) 



C'(0) = -ln\/2^. 



D 



Lemma 3.5. For fc G N, 

(3.11) Afe(i) - (-l)'=-ii?,2i-^ ln2 - (1 - 2i-'=)fcC'(l - k). 



In particular, 
(3.12) 

(3.13) A2r. 



A,a) = --ln2, 



.l^ ( ^,„+i (l-2-^")(2n+l)!C(2n+l) 

-1(2) - l-i) KTKZvi;; 



2(27r)2 



n > 1. 



Proof. Differentiate C(z, i) = (2^ — 1) C(z) at z = 1 — A; and use the identity 
(3.14) 



c(i-fc)= ^'^r'^' . fceN. 



n 



Hurwitz's Fourier representation ( |1.3D for C(z, (7) in the range < g < 1 and 
negative z is absolutely convergent and thus can be used to directly obtain a Fourier 
representation for the function Ak{q) in the range < q < 1. 

Define 



(3.15) Ciz,q) = Y, 



cos{2TTnq) 



and S{z, q) = y 



sm(2Tmq) 



so that 



^, , v-^ Inn 0* ri/ N V^ ^11"^ ■ /^ \ 

—-C[z,q)^-} — -cos(27rnq) and --5(z,g) = -> — — sm(27mq). 

n— 1 n— 1 



dz 
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Now replace z hy 1 — z in the Fourier representation (1.2) and differentiate to 
produce 



2r(z) 

(3.16) C'(l-.,?) = ^x 



■ , , TTZ TT nZ 

*(z)cos— + -siny 



C{z, q) + ^*(z) sm Y - 2 ^°^ Y 



S{z,q) 



TTZ U TTZ U 

COS -:— T^C(z, g) — sin -—-—S{z, q) 



2 dz 



2 92 



where ^(z) := ln27r — ^p{z). 

For z a positive integer, the f unct ions S (z, q) and C(z, q) are related to the 
Bernoulli polynomials (in view of (1^) and (L6)) and the Clausen functions. The 
latter are defined by 



(3.17) 

and 
(3.18) 

One has 



Chnix) = J2 



sin kx 



-, n>l 



fe=i 



oo , 

CWi(.^) = E^^, n>0. 



fe=i 



j^2n+l 



n— 1 ^ "^ 

(3.20) C(2m + l,g) = f^ £^g^ ^ CT2™+i(27rg), 

n=l 

(3.21) 5(2m + 2,g) = E ''"^,l"+7^ - g^2.n+2(2^g), 

n=l 

,,,,, r^^o ^o ^ ^ cos(2^ng) (-1)" (27r)^"+^ ^ 

(3.22) C(2m + 2,g) = E ^2^+2 = 2(2m + 2)! ^^-+^(g)- 

The value z — 2m + 1 yields, upon using '>Jj{k + 1) = —7 + Hk, the expression 

(3.23) A2rn+i{q) = {H2m ~1~ In 2^) B2™+i(g) 



+ (-1)' 



,2(2m+l)! 



(27r)2™+i 
Similarly, for z = 2m + 2 we find 

(3.24) A2,n+2{q) = (H2rn+1 - 7 - In 27r) B2m+2{q) 

+ (-1) 



Elnn . ,„ , TT v-^ cos(2TTnq) 
^-—p sin 27mq) + - > , ^, ' 

,n— 1 n— 1 



+i2(2m + 2)! 
(27r)2'"+2 



Elnn \ "' v^ sin(27rnq) 

-^^^ C08{2TTnq) - - }_^ „2m+2 



2 ^^ n- 

11=1 
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Lemma 3.6. For Q < q < I the function ^2(9) is given by 

y-ln2^)(g2_5+i) 

Emn , . 1 v-^ 

-^ C0si2nnq) + t^^^ 



(3.25) A^iq) = (l-7-ln2^)(g2_5+i 



n— 1 n— 1 



In particular, 



^2(1) = 2C'(-1), 

A2{\) - -C'(-l)-j^ln2, 

^2(i) = -K'(-i) + ^, 

where G is Catalan's constant. 



Proof. The expression ( 3.25 ) foUows directly from ( 3.2-^ ) at m = 0. D 



Special values of the derivative of the Hurwitz zeta function at the rational 
arguments g = ^, I, 4, |, i and |, for fc odd, have been given in [ni|. Special values 



of A};{q) for q > 1 can be obtained by using Lemma 3.3 a sufficient number of times. 

Example 3.7. 

^2(2) = 2C'(-1), 

^2(3) == 2C'(-l) + 41n2, 

^2(1) = -C'(-l)-iln2, 



^2(1) = -iC'(-l) + --ln2. 



G_ 

2^ 



Several of the indefinite integrals derived in the next sections can be conveniently 
expressed in terms of a family of functions closely related to the negapoly gamma 
family introduced by Gosper [p| as 

V'-i((z):=lnr(g), 

(3.26) /•« 

il^-k{q) := / tp-k+i{t)dt, k>2. 
Jo 

These functions were later reconsidered by Adamchik |l| in the form 

1 ri 



(3.27) ^-'»('^) = ^fc3^y (g-t)'-2lnr(t)di, k>2. 
We introduce the balanced negapoly gam,m,a functions by 

(3.28) V^'-'-Hg) := Y, [Ml) - Hk^iBkiq)] , 

where fc G N and Hr is the harmonic number [Hq :— 0). For instance, 



(3.29) V(-i)(g) = A,{q)=\nT{q) + C{0), 

(3.30) ^^-'Hq) = lA,{q)-^B,iq). 



10 
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Note. We can express the derivative of the Hurwitz zeta function at non-positive 
integers as 

(3.31) C'(_^,5)^H^(-i-'-)(9) + -^i3i+,(g), r e Nq. 

1 + r 



Lemma 3.8. The balanced negapolygammas can be expressed as 



(3.32) 



^(-'=)(g)=e-''^^ 
oz 



=7Z 



az,q) 



r(i-z) 



= l-/c 



Proof. Perform the derivative and use V'(fc) = Hk~i —7, in addition to the definition 
(O) and the identity (O). D 



We shah now study some of the properties of the balanced negapolygamma func- 
tions il^'^^^^iq). The adjective balanced is motivated by the following result: 



Lemma 3.9. For fc e N, 
(3.33) 



lp'^-^\q)dq = {). 



Proof. Use (3.6) and the analogous result for the Bernoulli polynomials, 

/ Bu{q)dq^O, 
Jo 



for k eN. 



D 



Lemma 3.10. For fc e N, 



(3.34) 



^^(-'=)(,)=^(-'=+l)(q). 



Proof. For fc > 2, (3.34) is a direct consequence of Lemma 3.4 and the well-known 
property of the Bernoulli polynomials, 



(3.35) 



— Bkiq)^kBk-iiq). 



The case fc = 1 follows directly from ( ^.29 ) since 

d 



dq 



lnr(g)=^(<z)-V'"H'Z)- 



D 



The functions tp^~''\q) defined by (3.28) are thus closely related to the ones 
introduced by Gosper. In fact, the precise relationship between both families of 
polygammas of negative order is 

'=-1 fe-l-r 



(3.36) ^^-'^\q)=^jMq) + J2 



r!(fc-l-r)! 



[a-r)+HrC{-r)], 



where we have used the evaluation 

(3.37) ^(-i-'-)(0) = i[C'(-r)+i7,,C(-r)], 



eN. 
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Note. According to ||l[, 
(3.38) C'(-r) 



HrCi- 



In Ar 



where the Ar are the generahzed Glaisher constants defined by Bendersky ||^. 

Lemma 3.11. The balanced negapolygamma functions i/"^ '^•'(9) admit the follow- 
ing Fourier expansion in the range Q < q < 1: 



(3.39) ^(-''(9) 



{2nf 



v^ln(27rn)+7 fcTr 

2^ ^fe cos(2™q - — ) 



TT .^-^ 1 . kiT 

ol^Tik M^^riq - — ) 



2 ■'—' n 

n—l 



Proof. In the definition ( 3.2^ ) of 'tf/^^^q) substitute, depending on the p arity of fc , 
the Fourier expansions ( 3.23 ) or ( 3.24 ) for the functions Ak{q) and ( 3.19|) or ( 3.22 ) 
for the Bernoulli polynomials. D 



Lemma 3.12. For fc e N, 

(3.40) v^"'"' (g + 1) = V^^"''^ (q) + 



^k-l 



[Inq-Hk-i]. 



(fc-1)! 
Proof. In the definition ( ^.2^ ) use (^^) and the property 

B^{q+l) = B„M+mq"'-\ 
satisfied by the Bernoulli polynomials. 

Corollary 3.13. For r e N, 

(3.41) ij^-^-^Hl) = ^^-^-^\Q) = i [C'(-r) + HrCi-r)] . 

Proof Set fc = 1 + r and evaluate (|3.40| ) at q = 0. Then use (|3.37[) . 
4. Indefinite integrals of polygamma functions 



D 



D 



Integral formulae involving the polygamma functions can be obtained from the 
corresponding ones for ^(2,(7), like (2.3) or ( p.5[ ), by taking the limit z ^ m = 
2,3, ... . This limit is not trivial in general, because of the vanishing of the Pochham- 
mer symbol (1 — z)j^i for j > m — 2, and the appearance of the function C(l, q) 
on the right-hand side of the formulae cited above. Similarly, differentiation at 
2; = 1 — fc leads to evaluations involving the functions Ak (q) and the negapolygamma 
functions. Because of the connection ( |3.29| ) between ?A^^^-'(<z) and lnr(<7), and of 
the latter with the function lnsin(7rg), we also obtain indefinite integral formulae 
involving these last two functions. 

The next theorem gives the moments of the polygamma functions in terms of 
themselves and the balanced negapolygammas defined in Section 0. The proof of 
the theorem rests on the following result: 

Lemma 4.1. Let p e Nq. Then 



(4.1) 



p 

E 

r=0 



(-ir 



p- 



qP-^Br+ija + bq) 

br+l 



-' + {-ir 



Bp+i{a) 
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Proof. The basic identity 



P+i 



(4.2) 



Bp+i{x + y) = Y^ 



r=0 



P+1 



Br{x)y 



p+1 — r 



appears in p2[ 19:5:3. The identity of the lemma is obtained by replacing x by 
a + bq and y by ~bq in (4.2). D 



Theorem 4.2. Let n e No, rri e N. Then 

_(-l)^ 



(4.3) Jq-^l,(^)(a + bq)dq = n\J2^ 



~n-j„l,{m~j-l) 



{n - j) 



:TTg""^V' 



(a + 6<7). 



Proof. In view of relationship (1.10), we set z = m + e in (2.3) and take the limit 
e — > 0. For a general value of n there will be three kinds of terms to consider in the 
smii on the right-hand side of (2.3): 
a) For < J < TO — 2 no singularities will arise and we simply have 



am^j+e,q) ^ (-l)"+\ ,,^_,_n 



lim 



V'^ 



(9). 



+0 (-m - e)j+i 

b) For J = TO — 1 we use the results 

(-m - e)ra = (-l)"m! [l + H„,e + 0{e^)] 

and 

1 



ai+e,?) 



^iq)+Oie' 



where ip{q) — V' (?) is the usual digamma function, to obtain 



C(m ~j + e, q) 



i-m - e)j+i 



(-1)™+! 



1 



— +7/^(g)+iJ„ + 0(e) 



c) For i > m, say j = m + r with r > 0, we use 

(-TO - e)™+i+, = (-l)'"+iTO!r!e [l + (iJ„ - i/,)e + 0{e^)] 
and 

C(-r + e, g) = ({-r, q) + eC{~r, q) + ©(e^) 
to obtain 



C(to -j + e, q) 



(-TO - e)j+i 



(-1) 



m+l 



j— m.+r 



7n! 



lBr+i(g) , „ Br+i{q) 



e (r + 1)! 



(r + 1)! 
-^(-i-'^)(g) + 0(e) 



where we have used the definition (3.28) of the balanced egapolygamma function 
and (1.6). When all the terms are added up, we find, in view of Lemma 4.1, that 
the coefficients of the 1/e singularity and the term proportional to the harmonic 
number Hm reduce to a g-independent constant, which can be dropped. This proves 
the theorem. D 



A similar result holds for the moments of the digamma function: 
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(4.4) 



3=0 



{-ly 



n-i,/,(-i-i) 



6J■+l(?^-j)! 



-tX^'^^ 



(a + bq). 



Proof. Use (1.9) in (2.3) and proceed along the same lines as the proof above. D 



Note: Result (4.4) can be identified as an extension of Theorem 4.2 to the case 
m — 0. 

Note: Adamchik m H| has provided the alternative representations 



r x-^pix) dx = i-ir [^f^ c'(-")) 

n 



fc=0 



:)="-' G"(-.=)-^^^i4^ 



and 



(-fc). 



Bk+i(z)Hk 
fc + 1 



f x-i^{x)dx^Y.^-lf'-^{t){^' 

Jo f.^Q \ / \ 

' E(-l)'fc! {l} l°gG.+r(. + 1) + (-l)"H„^^^±i^^pi^, 

where \ if are the Stirling numbers of the second kind and Gk{z) is the multiple 

Barnes function. The first representation can be directly obtained from (4.4) at a = 
and & = 1, by explicit evaluation of the integral between q = Q and q = z, making 
use of ( 3.37 ) for ^(^■'"^^(O) and of the definition ( 3.2^ ) of the negapolygamma 
functions. 



By differentiating formula (2^) at z = 1 — m we can arrive at the following 
result for the moments of the balanced negapolygamma functions, which extends 
Theorem 4.2 to negative values of m: 



Theorem 4.4. Let n e Nn. rri e N. Then 



(4.5) fq^^^(-"^){a + bq)dq^nlJ2T-ri 

3=0 



i-iy 



."-i,/,(— 'n-i-l) 



{n - J) 



:tt«""'V' 



{a + bq). 



Proof. Differentiate (2^) at z = 1 — to and use the result 
d 



dz 



(1 — z)j+i — — (to)j + i [Hm+j — H-m-l] 

z — l — m 

to obtain the following result for the moments of the function Am (q) : 



(4.6) 



q"A„,{a + bq) dq = m! n! ^ — ^ 



(-IP 



3=0 



6J+i(n- j)!(to + j + 1)! 



{H, 



m-\-j 



H 



7n — 1 



Am+j+iia + bq) 



)Bjn+j+iia + bq) 
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This result is equivalent to the statement of the theorem in view of definition (3.28) 
of the balanced negapolygammas and result ( ^.7[ ) for the moments of the Bernoulli 
polynomials. D 

As a particular case of the last theorem we obtain a formula for the moments of 
the loggamma function Inr(gf). 

Example 4.5. Let n G Nq and a, b eR. Then 

g"lnr(a + bg)dg = lnV2^^+7i!V ~ / ^ V^(-^--')(a + bg). 

n + 1 ■^ oJ+^(n — jj! 

This generalizes Gosper's result S, which establishes that all integrals of the form 
/ q" Ing! dq are expressible in terms of C,'{—i^ q), with 1 < j < n + 1. 

As particular cases of (^) we have 

g" In r(g) dq = In x/2^^ + n! ^ S^lA^^i-^-^ (^) 

and 

(7"lnr(l -q)dq = ln\/2^-^ - n\^-^——ijj'^-'^-^\l - q). 

n + ,^^ (n — J). 



Proof. Apply theorem (4.4) to V'*- ^Kq) = Inr(g) — In V2Tr. The expressions (4 



and (4.9) correspond to a = 0, 6 = 1 and a = 1, 6 = — 1 respectively. D 



The two special cases of Example |4.5| are now combined with the reflection 
formula for the gamma function 

(4.10) r(g)r(i-g) = -^ 

sm nq 
to obtain an expression for the moments of Insinvrg. 

Example 4.6. Let n G Nq. Then 

(7"+i In 2 



(4.11) Jq" 



In sin nq dq 



n+1 



Proof. Use the reflection formula for T{q), results ( 4~8| , 4~9| ) and the definition (3.28) 



to produce (4.11). The term that corresponds to the Bernoulli polynomials disap- 
pears in view of 

(4.12) {-iyB,+2{q) = B,+2(l-g). 

D 
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Example 4.7. 

(4.13) / e'^lnsinvrg dq = -e' 



ln2 + ^ 



{-iyA,+2iq)-A,+,il-q) 



(j+2)! 



Proof. Divide (4.11) by n\ and sum over n. 



D 



Example 4.8. Integrating (4.13) by parts yields 
(4.14) 



e"^ cotgTrg dq — 



In sin Trq + In 2 + ^ 

3=0 



{-iyA,+2{q)^A,+2{l-q) 
a + 2)! 



5. Some definite integrals 

Some of the definite integral formulae given in 0|, in the range (0, 1), can be 
obtained directly from the indefinite integral formulae given in Sections and m 

Example 5.1. Evaluate equation (^.3[) between and 1 to obtain 



(5.1) / q-a^,a + bq)dq^n\J2 ' >^' / 



i-iyC{z-J- l,a + b) 

(n-jy. 



n\{-iy 



{C{z — n — 1, a + fe) — C(z — n — 1, a)) . 



6"+i(l-z)„+i 
As a particular case we obtain formula (12.2) of H: 

Corollary 5.2. Let n e No and z eR, with z -n -1 <0. Then 



(5.2) 



q C[z, q) dq^n\y ^ --. 

^ (l-z)j+i(n-j)! 



Proof. Set 6 = 1 in (5.1) and use the identity (3.4) and the hypothesis z — n— 1 < 
to get rid of the last term in (5.1) in the limit a — > 0. D 



The evaluation of formula ( 4.11 ) between q = and q = 1 leads to formula (5.6) 



of0 

Example 5.3. Let n e Nq. Then 



(5.3) 



LfJ 



'E' 



(-l)'=C(2fc + l) 



' n, , ■ N . In2 

a [msiiiTT q) da = \- nl , , ^^, , 

^ ^ ^' ^ 71+1 ^^ (27r)2fc (n + 1 - 2k)\ 



Proof. Direct evaluation of the right-hand side of (4.11) gives, in view of property 



(5.4) 



g" In sin nq dq — 



In 2 

n+l 



71 — 1 ^ 



Clearly only the terms with j odd, say j = 2fc — 1, survive in the sum. (5.3) then 
follows directly from (3.9). D 
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Example 5.4. Let n G Nq. Then 



.1/2 

(5.5) / q"-\ii{smTrq)dq ^ 
Jo 



1 



2»+i 



In 2 

n + 1 



E 

fe=i 



(-l)''(22^-l)C(2fc+l) 
(27r)2'=(n + l-2/s)! 



1 - f-D" 



Proof. Use ( |3.1lD in ( |4.11 
For instance, 



D 



1/2 



(7ln(sin7r(7)dg 





1/2 



q ln(sin7r(7)(ig 



8 ^ 16^2' 

-lln2 + ^, 
24 167r2 ' 



ri/2 1 

/ (7^1n(sin7rq)dg = -—In 2 
Jo 64 



9C(3) 93C(5) 
647r2 1287r4 ' 



Example 5.5. Formulae (|4.7| ) or ( [4.8| ) allow us to derive Gosper's formulae for 
integrals of Inr(g) |g] in a very economical way. For instance, setting n — 0, 
a = b= I'm (O) yields 



(5.6) J' lliTiq + l)dq = q\liV2^+ ^A2{q + 1) - ^^2(9 + 1) - C'(-l) + ^^2 

Evaluation of the right-hand side at q = i and j yields, respectively, 

/■1/2 

/ \nT{q+l)dq = 
Jo 



8 24 2 2^ ^ ^' 



\nT{q + l)dq = - A _ ^ 1^2 + 1 In V2^ - ^C'(-l) ^ 







32 2 



47r' 



which can easily be seen to be equivalent to Gosper's formulae, after using Rie- 
mann's functional equation for the Riemann zeta function to express C'(— 1) in the 
form 



(5.7) 



/■'(o) 1 , 

C'(-l) = ^-TTT(21nV2;^ + 7-l) 



27r2 12 



Example 5.6. For fc, k/ G 



tp^-''\q)'ip^~^'\q)dq 



2cos(/c — k'' 

(27r)'=+fe' 



C"(fc + k') - 2(7 + ln27r)C'(fc + k') 



+ (7 + ln2^)2 + _ C(fc + fc') 
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The special case k = k' ~ 1 reduces to 



7 V ^V 1 , A^ 4 , 2 



(InTiqjydq = ±^ + — + -j\nV2n + -\nW2n 



I ^o^ /^,C'(2)^C"(2) 



given m 



Proof. Use the representation ( 3.32 ) of the balanced negapolygammas to obtain 
the desired result, by direct differentiation of the formula 

(5.8) j^ C{z\ q)az, q)dq = ^^^^LL' C(2 - z - z') cos [ ^ ) - 
valid for real z, z' < 0, given in 0. D 
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